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Abstract. We study the problem of two superconducting quantum qubits coupled via a 
resonator. If only one quanta is present in the system and the number of photons in the 
resonator is measured with a null result, the qubits end up in an entangled Bell state. Here we 
look at one source of errors in this quantum nondemolition scheme due to the presence of more 
than one quanta in the resonator, previous to the measurement. By analyzing the structure 
of the conditional Hamiltonian with arbitrary number of quanta, we show that the scheme is 
remarkably robust against these type of errors. 
PACS:03.67.Lx,85.25.Cp,74.50.+r 



Superconducting qubits have emerged in the last decade as a reliable benchtop onto 
which quantum-information experiments can be implemented. Architectures that demonstrate 
entanglement at the level of two qubits have been proposed and demonstrated experimentally 
(for extensive references, see pQ). Here we concentrate on a scheme that has been proposed in 
[T], namely entanglement of the two qubits via nondestructive measurements of the number of 
quanta in third quantum system coupled to the two qubits. 

This third quantum system can be for example a third qubit, a Josephson junction, an 
electromagnetic resonator, or a nanomechanical resonator. Schemes such as these have been 
proposed initially in quantum optics p], and, together with interaction-free experiments [3] 
and partial-collapse measurements [3], they have no classical analogue as they explore the 
consequences of the projective measurements in quantum mechanics. If the number of excitations 
present in the system is 1, there is no principial difference between a two- level system and a 
resonator as the intermediating system. However, in the case of a resonator as a connecting 
system, it can well be that the number of excitations present in the system, n, is larger than 1. 
We analyze here precisely this situation. 

The Hamiltonian of two qubits (biased at the optimal points) coupled via a resonator is 
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In the following, we take the two qubits identical and on-resonance with the coplanar waveguide 
resonator, Ej\ = Ej2 = u r , and the same values for the qubit-resonator coupling constants 
91= 92= 9- 

We now consider the situation in which the resonator is out-coupled to a detector. As 
a detector, one can use a single large Josephson junction and perform switching-current 



eigenvalue 


eigenvector 










(0,-1,1, Of 






-ig^2{2n - 1) 


1 n-l 1 In - / 2n-l „■ / 2n-l 
V2(2n-1) ^ V"- 1 ' ^V 2 ^- 1 )' V 2 (»-l)' J 






ig^/2{2n - 1) 


1 n-l 1 1 n ■ 1 2n-l „• / 2n-l i\ T 
V 2(2n-l) I, V"-!' V 2 (n-l) ' 1 V 2(n-l) ' ^ 



Table 1. Eigenvectors and eigenvalues of —iH c for k = T = 0. 



experiments [5J. Under the condition that no event is detected by the junction, and including the 
dissipative terms ri = T2 = T for the qubits and the resonator, the conditional non-Hermitian 
Hamiltonian [6], in the interaction picture, reads: 
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We notice that the conditional Hilbert space of the system splits naturally into 4-dimensional 
subspaces with n excitations distributed between the two qubits and the resonator. The 



Hamiltonian H c projected onto the subspace spanned by the basis states |1. 
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This non-Hermitian matrix, —iHc , for fc = T = 0, has two eigenvectors with eigenvalue 0, 
and two other with nonzero eigenvalue (see Table [T]) . We evolve now an initial state with n 
excitations |^(0)} by the conditional Hamiltonian and find the state at any time t under 
the condition that the resonator did not decay, exp(— iHc t)\ij)(0)). The results are shown in 
Fig. Q. 

From Fig. ([1]) we see that indeed for n = 1 the probabilities corresponding to the states 
|2,n = 1) and |3, n = 1) are equal. It is in fact known [TJ [2], that, if the system starts in 
the states |3,n = 1) or \2,n = 1} in the asymptotic state the qubits will be projected to a 
maximally entangled Bell state, |3,n = 1) — |2,n = 1)+ = (1/a/2)|0) p (g> (| fl) — I IT)) (we have 
verified the also the minus sign numerically). This corresponds to the second eigenvalue of the 
nondissipative Hamiltonian from Table [TJ 

In the case of n / 1 photons a completely different state results, namely the asymptotic 
state corresponds to the first eigenvalue in Table UJ This statement is valid no matter which of 
the four states (or combinations of them) |l,n), |2, n), |3, n), |4, n) is taken as the initial one, 
although of course the intermediate-time evolution is different. 

From Fig. ([1]) we see that the ratios of the two surviving probabilities are indeed 
fij n \/Pa n \ = (n — l)/n; we have also checked numerically the signs. The asymptotic state 
is then 
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Figure 1. The conditional occupation probabilities for the four states |1, n), |2,n), |3,n), and 
|l,n), for n = 1,2,3,4 and initial state |2, n). 



Interestingly, in the limit of large n in which we can take \n) p ~ |n — 2) p and factor out the 
resonator, this state results in another two-qubit Bell state, (| jj) + | ||})/\/2. It is intuitively 
clear that due to the presence of extra excitations in the system the probability of decay will 
increase and there will be no asymptotic stationary state as in the case n = 1. To understand 
this quantitatively, we calculate the probability to have no photon emission until time t. This is 
given by the norm of the wavefunction evolved with the non-Hermitiam conditional Hamiltonian 
with n excitations, 
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This probability is plotted in Fig. [2]for n = 1, 2, 3, 4. To make the interpretation straigthforward, 
we took ri = r2 = r = 0, in other words the only possibility of decay is through the resonator. 
We see that for n = 1 and an initial state |2,n = 1) (or equivalently |3,n = 1)) the probability 
Pq U ~ 1 ^ goes asymptotically to 0.5, as it should in the case k 3> T, as was found in [U [2], 
corresponding to an eigenvalue of the form (| j j) + | j j)) /\/2. For the other possible initial state 
with n = 1, namely |l,n = 1), the probability decays fast, therefore the errors introduced by 
these states will be asymptotically suppressed. 

For n / 1 the probability of no decay goes very fast to zero: this means that the efficiency 
of the scheme is much smaller than in the case n = 1, in the sense that there will be much 
less "favorable" events in which no decay of the resonator is detected. To make things more 
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Figure 2. The probability of no photon emission in the limit T 3> k and for different 
number of excitations. The graph on the left corresponds to an initial state |2, n = 1), while 
that on the right to \l,n = 1). The other parameters were g\ = g% = 100 MHz, k = 50 MHz, 

ri = r 2 = o . 



formal, suppose we start with an initial mixture of states |2, n = 1) (desirable, probability P^) 
and thermal excitations or other kind of errors and fluctuations giving n/1 excitations in the 
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system with probabilities Pex ■ The explicit form of this initial state does not matter, because 
the asymptotic result is the same: the asymptotic state will be a mixture of the Bell state we 
want to obtain, (| ][) + | []))/\/2, which comes with probability 1/2 and states corresponding 

to \ipasym), which come with (the much smaller) probabilities Pq 1 ^, 
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In conclusion, irrespective to the initial weights P^\ Pex\ this density matrix approaches 
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exponentially fast the desirable Bell state, due to the fast decay of the probabilities Pq ■ Of 

course, the initial state and p(°\ P^ do have a role to play in the efficiency of this process 
(how many favorable non-counting events we get), but not in the structure of the final state 
under the condition of no de-excitation events detected. 
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